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Abstract: Neural tube defects (NTDs) are congenital anomalies resulting from the failure of neural 
tube closure during embryogenesis. The precise molecular mechanisms underlying this 
multifactorial disease is poorly understood, although single nucleotide polymorphisms in genes 
involved in the one-carbon metabolism cycle are believed to contribute towards NTD 
development. Among them is 5-methyltetrahydrofolate-homocysteine methyltransferase 
reductase (MTRR). Protein function prediction algorithms (PolyPhen-2, PROVEAN, SIFT, SMART- 
Ensembl) were employed to evaluate its pathogenicity potential caused by the replacement of 
isoleucine with methionine. Seven NTD patients and 12 of their parents were recruited for this 
study. DNA samples were collected through blood or saliva whereby the extracted DNAs were 
then sent for whole exome sequencing (WES). Zygosity of the variant was confirmed from WES 
data of each subject and further validated through polymerase chain reaction (PCR) and Sanger 
sequencing. The results revealed that 57% of patients and 83% of parents carried rsl801394 
mutation in their MTRR gene, based on either homozygous (G/G) or heterozygous (A/G) 
genotypes. Bioinformatics analysis of this missense mutation predicted that this change is 
damaging to MTRR protein function by 2 of the 3 predictor algorithms and that the change from 
isoleucine to methionine amino acid affects flavodoxin domain of the protein. This impacts 
enzyme activity within the one-carbon metabolism pathway, which is linked to the aetiology of 
NTDs. From population databases, this variant was considered common with a MAF >0.3, 
however, it was not found in the Singapore Genome Variation Project (SGVP), whose population is 
a closer representation of the Malaysian subjects investigated here. Hence, we explored the 
prevalence of this variant in other studies and found that its association with NTDs differed across 
populations worldwide. Finally, we conclude that rsl801394 may be an NTD risk factor in the 
Malaysian population and should be further investigated as a potential prenatal screening tool. 
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1.0 INTRODUCTION 

Neural tube defects (NTD), encompassing both brain 
and spinal cord, are congenital anomalies affecting on 
the average, 0.5 to 10 per 1000 pregnancies globally 
( Greene & Copp, 2014 ). NTDs result from failure of 
neural tube closure during embryogenesis, with the 
exposed neuroepithelial cells susceptible to 
degeneration in utero ( Mohd-Zin et al., 2017) . 
Depending on the site of the lesion along the neuraxis, 
the severity and types of NTD vary, with defects at the 
brain (anencephaly) being the most severe phenotype 
and incompatible with postnatal life. As for lesions 
affecting the spine, also known as spina bifida, lifelong 
neurological handicap in patients is to be expected. 

NTD is a complex disease with multifactorial aetiology, 
involving both genetic and environmental 
contributions. Research has been focused on 
elucidating the mechanisms and causes of neural tube 
closure failure, with periconceptional folic acid 
supplementation singled out as a key nutrient 
recommended by World Health Organization for 
primary prevention of NTD. As such, NTD cases in 
countries where fortification of folic acid is mandatory 
reports of reduction in prevalence ( Kondo et al., 2017) . 
While the recommendation of folic acid 
supplementation for women of child-bearing age 
periconceptionally has decreased NTD incidence by 
70%, the remaining 30% remained elusive ( Blencowe 
et al., 2010 ). For this 30%, it is postulated that 
predisposing genotypes work together with non- 
genetic factors in modulating neural tube closure 
( Copp et al., 2015 ). 

In view of folic acid's role in successfully mitigating a 
large proportion of NTD cases, it logically follows that 
genes related to folate metabolism are most widely 
investigated. The most promising candidate gene to 
date is 5,10-methylenetetrahydrofolate reductase 
(MTHFR) gene. It is known that MTHFR influences 
folate status by irreversibly converting 5,10- 
methylenetetrahydrofolate (5,10-methylene THF) to 5- 
methyltetrahydrofolate (5-methyl THF) and that the 
common C677T gene variant reduces its enzymatic 
activity. As a result, carriers of the T allele(s), both 
homozygous and heterozygous, have elevated 
homocysteine levels and lowered folate 
concentrations, contributing towards NTD 


susceptibility ( Kirke et al., 2004 ). The second 
polymorphism of this gene, A1298C, is suggested to 
decrease MTHFR activity through gene-gene 
interaction with C677T variant ( Zhang et al., 2013) . 
However, despite the strong case for C677T as an 
independent risk factor for NTD, it appears that the 
prevalence of these MTHFR variants vary across 
populations and ethnic groups ( Zhang et al., 2013; van 
der Put et al., 1997; Amorim et al., 2007; Del Gobbo et 
al., 2018) , including a negative association in Malaysian 
Malays ( Hayati et al., 2008 ). 

As such, research to probe for single nucleotide 
polymorphisms (SNP) in other genes in the one-carbon 
metabolism cycle is currently underway. One of the 
prime candidates is the 5-methyltetrahydrofolate- 
homocysteine methyltransferase reductase (MTRR) 
gene. MTRR gene, encoding for methionine synthase 
reductase protein, is targeted in light of its chaperone 
role for 5-methyltetrahydrofolate-homocysteine 
methyltransferase [MTR) protein, and the key 
regulatory role it plays in folate metabolism. Located 
on chromosome 5, a missense mutation in MTRR is 
believed to impair the functions of methionine 
synthase, an enzyme linking together folate and 
methylation cycles in the one-carbon metabolism 
pathway. A66G variant (rsl801394) has been reported 
to affect folate and vitamin B12 status ( Hiraoka & 
Kagawa, 2017 ), both of which are essential substrate 
and cofactor for the successful functioning of the one- 
carbon metabolism pathway. Based on the interrelated 
and interconnected reactions, perturbations in this 
network during embryogenesis may induce serious 
detrimental effects, including NTD ( De Castro et al., 
2010 ). A mouse study has also revealed that 
hypomorphic mutations in the MTRR gene of either 
maternal grandparent resulted in NTDs in their grand 
progeny, regardless of the MTRR genotype in the 
affected litters, and that this effect was transmitted up 
to 4 generations ( Padmanabhan et al., 2013) . This 
finding implies that a reduction in MTRR gene function 
has devastating consequences for subsequent 
generations; hence parental MTRR genotypes are 
worthwhile to be considered as well. 

The objective of this study, if a correlation between 
MTRR SNP and NTDs is established, may be suggested 
as a prenatal screening tool for the prevention of 
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NTDs. In addition, if future studies can demonstrate 
that adequate or high folic acid consumption is able to 
overcome the effects of MTRR A66G variant, the 
results from this study may serve as a basis for 
mandatory folic acid fortification in staple foods in 
Malaysia. 

2.0 MATERIALS AND METHODS 

2.1 Patient selection / Study subjects 

24 Malaysians were recruited for this study. From here, 
the 7 patients with sporadic NTDs were grouped into 
distinct families, comprising of 5 complete trios 
(mother-father-patient) and 2 single-parent families 
(mother-patient). There were no descriptions of NTD 
cases in other family members and none of the 
patients were the result of consanguineous parents. 
The patients were presented with spina bifida with 
neurological deficits, encompassing both spina bifida 
occulta and aperta. The remaining 5 subjects were 
healthy individuals unrelated to the patients and were 
included as normal controls. Blood or saliva samples 
were collected from participants and all provided 
written consent and ethical approval was given by the 
Medical Ethics Committee (MEC) in University of 
Malaya Medical Centre (UMMC; MEC reference 
number: 914.5). 

2.2 DNA extraction 

DNA was extracted from either peripheral whole blood 
samples or saliva of subjects. 5 mL of the venous blood 
sample was collected and lysed in a 5X red cell lysis 
buffer followed by the addition of proteinase K in SNET 
buffer. For saliva collection, Oragene DNA OG-500 and 
OG-575 kits (Genotek, Canada) were employed. The 
collection tube was incubated at 50°C overnight and 
subsequent DNA extraction steps were performed 
according to the manufacturer's protocol. 
Concentration and purity of extracted DNA were 
assessed with NanoDrop 8000 (Thermo Scientific, USA) 
spectrophotometer. 

2.3 Whole exome sequencing and computational 
analysis 

Five micrograms of each DNA sample was sent to Axeq 
Technologies for whole exome sequencing (WES). DNA 
library from the sample was prepared with lllumina 
TruSeq kit then processed for paired-end sequencing 
on the lllumina HiSeq 2000 system. Through the 
integration of tools, the reads were mapped to 
Genome Reference Consortium Human Build 37 
(GRCh37), also known as hgl9. MTRR rsl801394 
variant was confirmed from each subjects' exome data. 


Minor allele frequency was obtained from 4 population 
databases, namely Genome Aggregation Database 
(GnomAD), Trans-Omics for Precision Medicine 
(TOPMED), 1000 Genomes (1000G) and Singapore 
Genome Variation Project (SGVP). Next, the 
pathogenicity of the variant was analysed with the 
following protein function prediction algorithms- 
Sorting Intolerant From Tolerant (SIFT), Protein 
Variation Effect Analyzer (PROVEAN), Polymorphism 
Phenotyping v2 (PolyPhen-2) and Simple Modular 
Architecture Research Tool (SMART). 

2.4 Genotyping and result analysis 

The SNP data from WES were validated through 
polymerase chain reaction (PCR). PCR reaction was 
carried out in a final volume of 25 pL, comprising of 
lOx DreamTaq™ Green Buffer (Thermo Fisher 
Scientific®), 25 pM MgCI 2 , 20 mg/pL bovine serum 
albumin (BSA), 10 pM dNTP, 10 pM of each forward 
(5'-TGAAGTGATGAGGAGGTTTCTGT-3') and reverse (5'- 
TGTGGTGGTATTAGTGTCCTTT-3') primers, 50 ng 
genomic DNA, 5 pL DreamTaq DNA Polymerase and 
Milli-Q® ddH 2 0 (Merck, Germany). 

PCR reactions were performed in a thermal cycler with 
the following conditions: 4 minutes of initial 
denaturation at 94°C, followed by 35 cycles of 
denaturation at 95°C, the annealing temperature of 
59°C and elongation at 72°C. Each of these 3 steps was 
performed for 1 minute, followed by a final elongation 
at 72°C for 7 minutes. PCR products were resolved on a 
1.2% agarose gel and sent to First BASE Laboratories 
Sdn Bhd (Selangor, Malaysia) for Sanger sequencing. 

Allelic frequencies and the percentages of subjects 
carrying variant alleles were calculated based on the 
number of individuals with homozygous or 
heterozygous genotypes over the total number of 
individuals. 

3.0 RESULTS 

3.1 Details of MTRR SNP investigated in this study 

MTRR variant rsl801394 is located on chromosome 5, 
exon2, resulting in a non-synonymous protein change 
from isoleucine to methionine at position 66 amino 
acid (Table 1). The missense mutation is a result of a 
change of A to G nucleotide at position 7870973. High 
variant allele frequencies ranging from 0.36400 to 
0.46778 (36.40% to 46.78%) were detected among 
populations surveyed in GnomAD, TOPMED and 1000G 
databases whereas none were found in the SGVP 
database. 
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The missense mutation is predicted using 
computational tools to have pathogenic consequences 
onto MTRR protein function by 2 of the 3 protein 
algorithms used (Table 2), based on cut-off scores from 
each computational tool. 

3.1.1 SMART-Ensembl schematic 

Methionine synthase reductase is composed of 698 
amino acids, as represented by the horizontal grey line 
(Figure 1). Domains are presented in boxes, and the 
pink box shows low complexity region. The dashed red 
line indicates the position of rs!801394 variant 


occurred within the flavodoxin domain, which may 
impact protein function. 

3.2 Exome screening of MTRR and electropherogram 

Three genotypes were found in the 24 subjects: A/A 
(wild type), A/G (heterozygous variant) and G/G 
(homozygous variant). We found that 57% (4 out of 7) 
of the NTD patients and 83% (10 out of 12) of parents 
carried variant allele(s) at rsl801394 whereas all 5 
controls were wild type for MTRR variant rsl801394 
(Table 3), as confirmed from exome data and validated 
by PCR sequencing (Figure 2). 


Table 1. Position, effect, and minor allele frequencies of the MTRR variant rsl801394. 


Position 

Annotation 

Base 

Amino acid 


MAF 




change 

substitution 

GnomAD 

TOPMED 

1000G 

SGVP 

Chr5: 

7870973 

NM_002454.3 

C.A66G 

P.I22M 

0.46778 

0.43376 

0.36400 

Not found 


The position is based on GRCh37.pl3. MAF from GnomAD, TOPMED and 1000G databases were calculated from samples 
sizes of above 2500 each. For SGVP database, the MAF was divided into 3 major population groups in Southeast Asia, M= 
Malay (n=98); 1= Indian (n=95); C= Chinese (n=99). 


Table 2. Pathogenicity prediction of MTRR variant rsl801394. 


Prediction Software 

SIFT 

PROVEAN 

PolyPhen-2 

Pathogenicity 

(Score) 

Damaging 

(0.002) 

Neutral 

(-2.25) 

Probably damaging 
(1.000) 

Cut-off scores for each software 

:: SIFT = 0.050, PROVEAN = 

-2.50, PolyPhen-2 = 0.95 



Table 3. Frequency count of genotypes of rsl801394 in MTRR gene in 19 subjects. 


Subjects 


Genotype 


Number of subjects with 

AA 

AG 

GG 

variant ailele(s) 

Cases 

3 (43%) 

1 (14%) 

3 (43%) 

4 (57%) 

Parents 

2 (17%) 

8 (66%) 

2 (17%) 

10 (83%) 

Controls 

5 (100%) 

0 

0 

0 

AA = wild type; AG 

= heterozygous mutant; Het = 

homozygous mutant. 


Total = 14 


4.0 DISCUSSION 

In this study, we investigated the association of MTRR 
variant rsl801394 with neural tube defects. 14 out of 
the 19 Malaysian NTD patients and their parent(s) 
recruited in this study were confirmed to carry 
mutated MTRR rsl801394 allele(s), based on whole 
exome sequencing and subsequent validation by PCR. 
In each family, at least one member was detected as 
homozygous or heterozygous variant allele carriers. 


Inferring from GnomAD, TOPMED and 1000G 
databases, the occurrence of this SNP is common with 
a MAF of 0.36 and above, based on the accepted 
convention that any mutations recorded in more than 
1% or 0.01 of the general population are considered to 
be common ( Chen et al„ 2011) , and thus unlikely to 
bear clinical significance ( Kryukov et al., 2007 ). 
Flowever, this SNP was not found in the SGVP 
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Figure 1. Schematic diagram of MTRR protein domain as predicted by SMART-Ensembl using amino acid sequences 


Wild type 


Homozygous 


Heterozygous 


GCTCACiJATTTCTT 
120 130 


rs1801394 

(A66G) 


A/A 




Figure 2. Reverse sequencing data for rsl801394. SNP positions are indicated by arrows. 


database, which is interesting as the population group 
genotyped in SGVP, consisting of 292 Singaporean 
Malay, Chinese and Indian individuals are a closer 
representation of the Malaysian patients in our study 
in terms of genetic ancestries ( Teo et al., 2009 ). This 
finding is also consistent in our 5 healthy controls, 
whom comprises of the 3 main races in Malaysia. In 
the medical genetics field, differences in allele 
frequency profiles across populations are biologically 
meaningful ( Bomba et al., 2017) . Thus, by contrasting 
the prevalence of MTRR variant rsl801394 in more 
than half of the patients in our study against the rarity 
of this polymorphism in its demographic, it suggests 
that this genomic variation may be a contributing 
factor towards NTD development among Malaysians. 
Likewise, the high distribution of this MTRR variant in 
the unaffected parents hints at an underlying 
epigenetic inheritance pattern that eventually 
disrupted the neural tube closure process in their 
offspring. 

Protein function prediction tools provide valuable 
information on the pathogenicity likelihood of amino 
acid substitutions induced by the variant. Through 
evaluating the extent of sequence conservation at the 
missense site and structural features, the functional 
impact of the variant is then categorised into threat 
levels according to each algorithm's predefined 


threshold scores ( Raychaudhuri, 2011) . The adoption 
of 3 prediction tools increased the confidence level of 
the results and to account for false prediction biases in 
cases where the score is close to the cut-off value 
( Dakal et al., 2017 ). Although PROVEAN estimated that 
the rsl801394 variant is neutral to the MTRR protein, 
SIFT and PolyPhen-2 reported the variant as damaging 
with the score of 1.000 by the latter indicates 
maximum probability, thus is alarming. Hence, SMART- 
Ensembl was utilised to shed light on possible 
functional consequences by studying the affected 
protein domain. 

Methionine synthase reductase protein is expressed as 
a 698 amino acid long protein with 3 identified domain 
families- flavodoxin, FAD and NAD. The rsl801394 
variant resulting in I22M substitution affects the 
flavodoxin domain, which is a cofactor binding site. 
Biochemical analysis revealed that the methionine 
amino acid replacement affects interaction and affinity 
between the protein and its redox partner- methionine 
synthase ( Olteanu et al., 2002 ). Consequently, 
decreased reductive methylation and reactivation of 
methionine synthase results in reduced enzyme 
efficiency, which impacts the generation of S- 
adenosylmethionine (SAM) molecule required for the 
methylation process. It is to be noted, however, that 
since methionine synthase reactivation is only required 
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once in 2000 cycles, the reduction in SAM availability 
can be tolerated, as corroborated by a study that the 
demand for SAM during neurulation of embryogenesis 
is adequate ( Leung et al., 2017 ). Nonetheless, as 
methionine synthase is the crucial link connecting both 
folate and methyl cycles, a 3-4 fold compromise on its 
enzymatic activity may lead to other downstream 
implications in the overall one-carbon metabolism 
pathway, which has vital roles in embryogenesis ( Xu & 
Sinclair, 2015 ). Although the precise mechanisms have 
yet been determined, perturbations in this pathway 
have been linked to NTD development ( Li et al., 2016 ). 

Based on the probable pathogenicity of rsl801394 in 
mediating NTD risk and to overcome the small sample 
size limitation in this study, we explored available data 
to estimate the extent of association between this 
variant and NTD patients in other populations. The 
MTRR mutant genotype has been suggested to 
increase the risk of NTD development in Canadians 
( Wilson et al., 1999 ), Han Chinese in Tianjin province, 
China ( Fang et al., 2018; Wang et al., 2015 ), Polish 
( Pietrzyk et al., 2003) , Hispanics ( Zhu et al., 2003) and 
English ( Doudney et al., 2009) . Conversely, no 
associations were found in the Dutch ( Brouns et al., 
2008; van der Linden et al., 2006 ), Irish ( O'Leary et al., 
2005) , Chinese in Liaoning province (Zhang et al., 2019 ) 
and central China ( Liu et al., 2014) . Interestingly, one 
study reported that the rsl801394 variant protects 
against NTD and reduces development risk in their 
Northern United Kingdom cohort ( Relton et al., 2004) . 
As with the association analyses on MTHFR C677T 
where geographical and ethnic differences led to 
discordant results ( Copp et al., 2013 ), the conflicting 
reports on MTRR support the multifactorial nature of 
NTD. This is in line with the hypothesis that 
polymorphisms in a single one-carbon metabolism 
gene appear to have but a limited role in overall NTD 
risk determination and that the sum of gene-gene 
interaction confers a more pronounced pathological 
effect compared to individual mutant alone ( Relton et 
al., 2004; Greene et al., 2009) . 

This study has several limitations. Firstly, it is limited by 
the lack of detailed clinical data from each patient, 
which would otherwise help correlate genotype- 


phenotype relationship with the pathogenesis of this 
variant. Secondly, a genealogical analysis might 
uncover distant relationships among NTD patients in 
our cohort that would aid in understanding the 
transmission pattern of this variant. Thirdly, a larger 
dataset would provide a more convincing 
polymorphism association. 

5.0 CONCLUSIONS 

In summary, we report the association of MTRR variant 
rsl801394 in 57% of NTD patients and a high 83% 
prevalence of this variant in NTD parents in Malaysia. 
In light of the inconsistent prevalence of rsl801394 
variant across worldwide NTD populations, it is thus 
important to identify population groups whose genetic 
profiles may predispose them to increased NTD risks. 
Even though the sample size in this study is small, it 
provides initial evidence of the need to screen for this 
polymorphic variant in a larger Malaysian cohort whom 
may benefit from the prenatal screening of MTRR 
A66G mutants. Besides, while gaps in knowledge 
remain between MTRR gene polymorphisms and 
prevention of NTD through folic acid fortification, the 
findings of this study are hoped to serve as a basis for 
public health intervention in Malaysia in the future. 
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